Ising versus XY anisotropy in frustrated R 2 Ti 2 07 compounds as seen by polarized neutrons 
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We studied the field induced magnetic order in R^T^Ov pyrochlore compounds with either uniaxial (R=Ho, 
Tb) or planar (R=Er, Yb) anisotropy, by polarized neutron diffraction. The determination of the local suscepti- 
bility tensor provides a universal description of the field induced structures in the paramagnetic phase 
(2-270 K), whatever the field value (1-7 T) and direction. Comparison of the thermal variations of x\\ and^ ± 
with calculations using the rare earth crystal field shows that exchange and dipolar interactions must be taken 
into account. We determine the molecular field tensor in each case and show that it can be strongly anisotropic. 

PACS numbers: 71.27.+a, 75.25. +z, 61.05.fg 



The pyrochlore lattice of corner sharing tetrahedra offers 
the best model of geometrical magnetic frustration in three 
dimensions. In rare earth pyrochlore titanates R2T12O7, frus- 
tration arises from the subtle interplay of energy terms such 
as single ion anisotropy, exchange and dipolar interactions. 
Depending on the balance between these factors, one may ob- 
serve spin ice (R=Ho, Dy)±^ or spin liquid (R=Tb) 4 behav- 
iors, or complex magnetic orders stabilized by first order tran- 
sitions (R=Er, Yb, Gd) 5,6 . A major role in the selection of the 
ground state is played by the single ion anisotropy that arises 
from the trigonal crystal electric field (CEF) at the rare earth 
site. At low temperature, it may result in strong axial or pla- 
nar anisotropy, depending on the wave function of the ground 
state and of the close CEF excited states, if any. 

It is now widely accepted that Ho, Dy and Tb titanates 
show Ising-like behavior, the magnetic moments being con- 
strained to align along the local <111> axes. The use of 
the Ising model for the canonical spin ices H02T12O7 and 
Dy 2 Ti 2 7 in a large temperature range is justified by their ex- 
treme axial anisotropy and very high gap between the CEF 
ground state and excited states 7,8 . In Tb2Ti 2 7 , where this 
energy gap is much smaller 9 -*^, the Ising model holds only 
at low temperature; above 200 K, a Heisenberg-like behavior 
is recovered^ since many CEF levels become populated. The 
spin liquid ground state in Tb2Ti207 transforms into antiferro- 
magnetic (AF) order under applied pressure and/or magnetic 
field 12,13,14 . In zero field, none of these Ising titanates evi- 
dences long range magnetic ordering (LRO) down to the low- 
est temperature. 

The Er and Yb titanates, where the easy plane for the 
CEF ground state is the plane perpendicular to < 1 1 1 >W 5 , 
have been labelled XY-type pyrochlores, although they do not 
present extreme planar anisotropy. Er2Ti207 was suggested 
to realize a model XY-type antiferromagnet, for which the- 
ory predicts fluctuation-induced symmetry breaking, leading 
to magnetic ordering. It indeed orders antiferromagnetically 
at 1.2 K 5,16 , but can be driven into a quantum disordered state 
by application of a magnetic field 17,18 . In Yb2Ti207, a tran- 
sition detected at 0.25 K 19 was shown by 170 Yb Mossbauer 
spectroscopy and yuSR 6 to occur towards a short range corre- 
lated state. This transition is first order and consists in a sharp 
and strong decrease of the spin fluctuation frequency. How- 
ever, a time dependent ferromagnetic order was also reported 



in Yb2Ti207 , and the application of a moderate magnetic 
field was shown to drive it into a stable LRO state 21 . 

In the pyrochlore lattice, selection between Ising, Heisen- 
berg or XY models cannot be based, as usual, on the analysis 
of the macroscopic properties of a single crystalline sample 
in a magnetic field because of the presence of four anisotropy 
axes, namely the <1 1 1> axes. Then only an average over the 
four local axes can be measured by classical methods, and 
no direct information can be obtained about the anisotropy of 
the local magnetic susceptibility of the R ion, nor of the lo- 
cal exchange/dipolar tensor. There is thus a clear motivation 
to get a direct access to these quantities and to quantify the 
axial/planar anisotropy in the rare earth pyrochlores. 

Here we show that the site susceptibility approach 22 devel- 
oped for polarized neutron diffraction is an accurate tool for 
obtaining the local susceptibility tensor^. We present our re- 
sults in two Ising-type (H02T12O7, Tb2Ti207) and two planar- 
type (Yb2Ti207, Er2Ti2C>7) compounds in a large temperature 
range (2K-270K). Knowledge of the ^-tensor allows one to 
determine the moment values and orientations in the field in- 
duced paramagnetic state, whatever the field direction! 4 -. The 
measured thermal variations of the components of \ are quan- 
titatively compared with calculations using CEF parameters 
deduced independently from inelastic neutron scattering spec- 
tra. To obtain good agreement with the data, one needs to 
introduce an anisotropic molecular field tensor, which takes 
into account both the first neighbor exchange and the dipolar 
coupling. 

Single crystals of H02T12O7, Tb2Ti207, Er2Ti207 and 
Yb2Ti207 were grown by the floating-zone technique, using a 
mirror furnace. The crystals were characterized by zero field 
neutron diffraction at 100 K and 5K. The nuclear structure 
factors were used to refine positional parameters, occupancy 
factors, isotropic temperature factors and extinction parame- 
ters, within the space group Fd3m. 

Neutron diffraction measurements were performed at the 
ORPHEE reactor of the Laboratoire Leon Brillouin. Po- 
larized neutron flipping ratios were measured on the Super- 
6T2 spectrometer 23 using neutrons of incident wavelength 
X n - 1 .4 A and polarization Po=0.98, and on the 5C1 spectrom- 
eter (/l w =0.84 A, Po=0-91), provided by a supermirror bender 
and a Heusler polarizer, respectively. For each compound, 
100-200 flipping ratios were measured, at selected tempera- 
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tures in the range 2 - 270 K, in a magnetic field of 1 T applied 
parallel to the [110] direction. The minimum number of flip- 
ping ratios used in the refinements with two parameters was 
30. The programs CHILSQ 24 were used for the least squares 
refinements of the flipping ratios. 

Assuming a linear response, the model assigns a suscepti- 
bility tensor^, of rank 3x3, to each crystallographically inde- 
pendent sit o 22 i 25 . Its components Xu depend on the site sym- 
metry, and the magnetic moment M d =^11 induced on a R ion 
at the 16d site in the unit cell is not necessarily collinear with 
the field, since the off-diagonal elements of x can be non-zero. 
The tensor ^ has the same symmetry as the tensor u describ- 
ing the thermal motion of atoms. The site susceptibility can 
be conveniently visualized as a magnetization ellipsoid con- 
structed from the six Xij parameters in much the same way 
as thermal ellipsoids are constructed from the Uij parameters. 
For a given ellipsoid, the radius vector gives the magnitude 
and direction of the magnetic moment induced by a field of 
1 T rotating in space. For the space group Fd3m, the symme- 
try constraints imply that x has only two independent matrix 
elements, xn andxn- At the 16 d site, the principal axes of 
the magnetization ellipsoid are oriented along the four local 
<111> axes. Their lengths are given by: x\\ = Xn + ^Xn and 
X± ~ Xn - Xn- These susceptibility components are deter- 
mined for each temperature by refining the flipping ratios. 

The thermal evolutions ofx\\ and^ ± are shown in Figs. 1-4 
for the four compunds. In the insets are represented the sus- 
ceptibility ellipsoids at 5 K (lower) and 250 K (upper), whose 
axes have been scaled by temperature to compensate for the 
Curie decrease, together with the magnetic moments induced 
by the field applied along [110]. The four R sites of a tetrahe- 
dron are then separated into two classes: the so-called a- sites 
(the two on the left), where the angle between the field and 
the ternary axis is 35.3°, and the /3-sites (the two on the right), 
where the field is perpendicular to the ternary axis 26 . The data 
were interpreted using a single ion approach, where the trigo- 
nal symmetry CEF and Zeeman Hamiltonian is written as: 
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<HcEF=Y J B n°n-gjfB J.H. 



(1) 



In this expression, the relevant coefficients are B^ v B®, B\, B®, 
# 3 and B^, and the O™ operators are expressed in terms of 
spherical harmonic a 27 ! 28 , which has the advantage that the B™ 
are roughly independent of the rare earth in an isostructural 
series. In equation (Q]), gj is the Lande factor of the rare earth. 
In order to take into account exchange and dipole- dipole in- 
teractions, we introduced an anisotropic molecular field tensor 
{/iy, A ± } and performed a self-consistent calculation, assuming 
a field component of 1 T either parallel or perpendicular to the 
ternary axis. 

The data for Ho2Ti207 are shown in FigQ] At low temper- 
ature, the ellipsoid elongations x\\ increase so markedly that 
they degenerate into needles, while x± is ver Y small, lying 
at the limit of experimental precision. The induced magnetic 
structure is such that the a- site moments point towards the 
field direction, while remaining along the local ternary axis, 
and the /3-site moments are zero. At high temperature, on the 
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FIG. 1: (color on line). Ho 2 Ti 2 7 : Susceptibility components;^ 
and^f ± versus T. The following CEF parameters (in K) were used 7 : 
B° 2 = 791, B° 4 = 3188, B\ = 971, B° 6 = 1007, B\ = -725 and B 6 6 
= 1178. Dashed lines: CEF only calculation; solid lines: calculation 
including effective exchange: A\\=0.05 T/fi B , A ± =0. For^ ± , both lines 
superimpose. 



contrary, the anisotropy is weaker and all the moments reori- 
ent close to the field direction. These features, and in par- 
ticular the peculiar shape of x±(T), are explained by the CEF 
level scheme of Ho 3+ (/=8, gj=5/4 ) in H02T12O7: the ground 
state is very close to the extremely anisotropic non-Kramers 
doublet \ J = 8;/ z = ±8), which has a moment of 10/i fi along 
the ternary axis and x± = 0, and the first excited state lies at 
about 250 K above the ground state 7,8 . The lines in FigQ] were 
calculated using the CEF parameters given in Ref |7|. We find 
that the CEF only curves (dashed lines) yield a good agree- 
ment with the data. Due to the smallness ofxAT), no reliable 
value for the molecular field constant A ± can be obtained. For 
X\\, introduction of a small positive (ferromagnetic) value A\\ = 
0.05(1) T/yUfi improves the agreement with the data (see FigQ]). 
Since the saturated Ho magnetic moment along <lll>is// = 
lOyUfi, this yields an effective field H\\ = A\\fi =0.5 T. The ex- 
change and dipolar interactions in the ground spin-ice state 
of H02T12O7 are well known at low temperature 29 , amounting 
to a positive effective exchange constant J e ff ^ 1.8 K, where 
the dipolar contribution dominates. This is equivalent to a 
field H e ff 0.27 T along the ternary axis. An exact compar- 
ison with our result is however difficult because the magnetic 
structure induced by the field in Ho2Ti207 is quite different 
from the zero-field moment arrangement in a spin-ice, hence 
modifying the dipolar field at a Ho site. 

The CEF anisotropy of Tb 3+ (J=6, gj=3/2) in Tb 2 Ti 2 7 
is considerably smaller than that of Ho 3+ in H02T12O7 
(see Fig At room temperature, Tb2Ti207 is practically 
isotropic 11,14,30 due to the relatively large number of popu- 
lated CEF levels 9 ^. Decreasing the temperature, it evolves 
progressively into an Ising type with a ratio of ellipsoid axes 
of the order of 10. A CEF only calculation using the parame- 
ters from RefUOl accounts practically for the^y data, and lies 
somewhat above the x± data. Introducing an anisotropic A- 
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FIG. 2: (color on line). Tb 2 Ti 2 07: Susceptibility components \\\ an d 
X± versus T. The following CEF parameters (in K) were used 10 : 
B° = 712, B° A = 3400, B\ = 1200, 5° = 1130, B\ =-700 and ^ 
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Dashed lines: CEF only calculation; solid lines: calculation 
including effective exchange: A\\=-0.05T//u B , A ± =-\T//u B . 



FIG. 3: (color on line). Er 2 Ti 2 07: Susceptibility components x\\ an d 
X± versus T. The following CEF parameters (in K) were used (see 
text): B° 2 =616, B° 4 = 2850, ^=795, 



£°=858, 



B\ = -493 and B\ 



=980. Dashed lines: CEF only calculation; solid lines: calculation 
including effective exchange: Ay =-0.15 T/fi B and A ± =-0A5 T/fi B . 



tensor improves the agreement, especially for x±, with the 
AF values: A ]{ =-0.05(2) T/jn B and A ± =- 1.0(2) T/ji B . We have 
checked that this i-tensor is coherent with the isotropic effec- 
tive value A - -0.33 T/ju B derived in Ref. 10 from the analysis 
of the high temperature powder susceptibility. For^ ± (r), the 
self-consistent calculation does not converge below about 5 K, 
which is the "ordering" temperature of the model; this is why 
there are no calculated points for^ ± below 5 K in Figj2] Both 
components of the i-tensor, which corresponds to an effective 
(exchange + dipole) low temperature interaction, are found to 
be of AF type, in line with the position of Tb2Ti207 in the 
phase diagram of RefBU Our own previous estimation 10 ten- 
tatively located Tb2Ti207 in the spin ice side of this phase dia- 
gram, but it was based on the high temperature isotropic deter- 
mination of A and could not take into account the anisotropy 
of the /i-tensor found in the present work. 

As seen from Fig(3j the anisotropy of Er 3+ (7=15/2, gj 
=6/5) in Er2Ti2C>7 is rather weak, and a spherical shape for 
the magnetic ellipsoids is recovered above 100 K. At low tem- 
perature, x± > Xh m agreement with the planar magnetic or- 
dering below 1.2 K found in Ref, 5. The CEF parameters for 
Er2Ti2C>7, extrapolated from those in H02T12O7, were refined 
so as to reproduce the energies of the two lowest CEF excited 
doublets^, at 6.3 and 7.3 me V, as well as the transverse com- 
ponent of the g-tensor of the ground Kramers doublet, which 
must match the in-plane saturated moment. The CEF param- 
eters listed in the caption of Fig{3] meet these requirements, 
with g||=2.6 and g ± = 6.8 for the ground doublet, the latter 
value yielding m sat = ^g ± \i B - 3.4 ji B , close to the measured 
value 3ju B . The CEF only calculation of the susceptibilities 
shows a good first order agreement for^y(r), and introduc- 
tion of the effective exchange improves greatly the agreement 
fo*X±(T)- We find the i-tensor in Er2Ti207 is also anisotropic 
and of AF type: A ± = -0.45(5) T/ji B and A l{ = - 0.15(l)T/// 5 . 
The in-plane effective exchange is thus stronger than its com- 
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FIG. 4: (color on line). Yb 2 Ti 2 7 : Susceptibility components;^ and 
X± versus T. The following CEF parameters (in K) were used (see 
text): B° 2 =750, £°=2350, £3 = 1200, £°=700, £^=-1000 and B 6 6 = 
700. Dashed lines: CEF only calculations, solid lines: calculation 
including effective exchange: A\\=2.5T/fi B , A ± =- 0.05T/// 5 . 



ponent along the ternary axis, which reinforces the XY char- 
acter determined by the CEF anisotropy. 

FigH shows the temperature dependence of x\\ an d X± f° r 
Yb 3+ (7=7/2, g/=8/7) in the planar anisotropy pyrochlore 
Yb2Ti207. Its magnetic ellipsoids are disk-shaped, like in 
Er2Ti2C>7, but a strong anisotropy persists up to much higher 
temperature because the excited Kramers doublets lie much 
higher in energy (700-1000K) 15 . In order to calculate x\\(X) 
and^ ± (r), we used CEF parameters extrapolated from those 
in Ho2Ti 2 07. The chosen set of parameters must also re- 
produce the thermal variation of the 4f quadrupolar moment 
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Q*f(T) oc QJ 2 z - J(J + measured by the y - y Per- 
turbed Angular Correlations technique 15 . The B™ coefficients 
listed in the caption of Fig@] meet these requirements and 
yield a rather good agreement with the experimental data for 
X±(T), implying A± 0. For^|(r), a ferromagnetic value 
A\\ =2.5(5) T/fiB is needed to reproduce the data. We checked 
that the powder susceptibility calculated using this i-tensor 
agrees with experiment, i.e. that it leads to an almost van- 
ishing paramagnetic Curie temperature below 20 K 15 . The 
present set of parameters, which differs from the previous set 6 
roughly by interchanging the B® and B® values, is coherent 
with those in the R2Ti 2 7 family. It yields a somewhat larger 
g\\ value for the ground doublet (2.25 vs 1.80), which implies 
that the Yb moment of 1.15yUg in the "slow fluctuation" phase 
below 0.25 K 6 probably lies along the <1 1 1> axis. 



R 


Ho 


Tb 


Er 


Yb 


h 


0.05(5) 


-0.05(2) 


-0.15(1) 


2.5(5) 


A ± 




-1.0(2) 


-0.45(5) 


0.00(5) 



TABLE I: Values of the two components of the molecular field tensor 
(in T/fi B ) in the R 2 Ti 2 7 compounds. 

In conclusion, polarized neutron diffraction allowed us to 
determine the local susceptibility tensor in the R2Ti2C>7 series, 



inaccessible by macroscopic measurements in single crystals. 
Its temperature dependence in the paramagnetic phase cannot 
be entirely accounted for by the sole crystal field anisotropy: 
a molecular field tensor A must be introduced, which encom- 
passes exchange and dipolar interactions, and which is found 
to be strongly anisotropic in most compounds. The values of 
the A components are gathered in Table I. Their sign and mag- 
nitude agree with the onset of planar AF ordering (Er) and 
spin ice short range order (Ho). For Yb, we find an extremely 
anisotropic ferromagnetic A- tensor along < 1 1 1 > with a mod- 
erate planar CEF anisotropy. By contrast, for Tb, which has 
a strong CEF anisotropy along <111>, the A- tensor is of AF 
type and characterized by a sizeable planar anisotropy. The 
anisotropy of the i-tensor is likely to be due mainly to the 
dipole-dipole coupling in the large moment systems (Ho and 
Tb), but our results suggest that exchange itself is anisotropic 
in the small moment systems (Er and Yb), where dipolar in- 
teractions play a minor role. This could be an important ingre- 
dient for theory in Tb2Ti207 and Yb2Ti207, where the ground 
state is not yet fully understood. 
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